Physical Models of Asteroids

Matematicko-fyzikální fakult

Distribution of spin-axes longitudes and shape elongations of main-belt asteroids

Context: Large all-sky surveys provide us with a lot of photometric data that are sparse in time (typically few measurements per night) and can be potentially used for the determination of shapes and rotational states of asteroids. The method which is generally used to derive these parameters is the lightcurve inversion. However, for most asteroids their sparse data are not accurate enough to derive a unique model and the lightcurve inversion method is thus not very efficient. Aims: To fully utilize photometry sparse in time, we developed a new simplified model and applied it on the data from the Lowell photometric database. Our aim was to derive spin axis orientations and shape elongations of asteroids and to reconstruct distributions of these parameters for selected subpopulations to find if there are some differences. Methods: We model asteroids as geometrically scattering triaxial ellipsoids. Observed values of mean brightness and the dispersion of brightness are compared with computed values obtained from the parameters of the model -- ecliptical longitude $\lambda$ and latitude $\beta$ of the pole and the ratios $a/b$, $b/c$ of axes of the ellipsoid. These parameters are optimized to get the best agreement with the observation. Results: We found that the distribution of $\lambda$ for main-belt asteroids is not uniform (in agreement with findings of Bowell et al., 2014, M&PS, 49, 95) and is dependent on the inclination of orbit. Surprisingly, the non-uniformity of $\lambda$ distribution is larger for asteroids residing on low-inclination orbits. We also studied distributions of $a/b$ for several groups of asteroids and found that small asteroids ($D<25\,$km) are on average more elongated than large ones.Comment: 10 pages; Accepted for publication in A&

Reflectance Spectra of Seven Lunar Swirls Examined by Statistical Methods : A Space Weathering Study

Higher magnetic field in lunar swirls is believed to deflect majority of incoming charged particles away from the lunar-swirl surfaces. As a result, space weathering inside and outside swirls should be different. We wanted to evaluate these differences, therefore we have examined seven swirl areas on the Moon (four mare and three highland swirls). We applied the Modified Gaussian Model to statistical sets of the Moon Mineralogy Mapper spectra. Using Principal Component Analysis (PCA), we were able to distinguish the old (weathered) material from both the fresh crater and swirl materials. The swirls did not follow the same behavior as the fresh material, nor were they fully separable. Additionally, we could distinguish between the mare and highland swirls (mare/highland dichotomy) based on the PCA and histogram plots of the albedo and strength of the 1000-nm absorption band. The mare/highland dichotomy can partially be caused by different FeO content in maria and highlands, which points to the existence of a threshold value that changes the spectral evolution due to space weathering. Slope behavior seemed to be dependent on whether the swirl was on the near- or far-side of the Moon, likely due to shielding of lunar nearside by Earth's magnetotail. Our results thus favor the solar wind stand-off hypothesis in combination with the fine dust transport hypothesis and point to the fact that micrometeoroid impacts generally do not reproduce the same weathering trends as all the space weathering effects together.Peer reviewe

Combining asteroid models derived by lightcurve inversion with asteroidal occultation silhouettes

Asteroid sizes can be directly measured by observing occultations of stars by asteroids. When there are enough observations across the path of the shadow, the asteroid's projected silhouette can be reconstructed. Asteroid shape models derived from photometry by the lightcurve inversion method enable us to predict the orientation of an asteroid for the time of occultation. By scaling the shape model to fit the occultation chords, we can determine the asteroid size with a relative accuracy of typically ~ 10%. We combine shape and spin state models of 44 asteroids (14 of them are new or updated models) with the available occultation data to derive asteroid effective diameters. In many cases, occultations allow us to reject one of two possible pole solutions that were derived from photometry. We show that by combining results obtained from lightcurve inversion with occultation timings, we can obtain unique physical models of asteroids.Comment: 33 pages, 45 figures, 4 tables, accepted for publication in Icaru

Sub-surface alteration and related change in reflectance spectra of space-weathered materials

Context. Airless planetary bodies are studied mainly by remote sensing methods. Reflectance spectroscopy is often used to derive their compositions. One of the main complications for the interpretation of reflectance spectra is surface alteration by space weathering caused by irradiation by solar wind and micrometeoroid particles. Aims. We aim to evaluate the damage to the samples from H+ and laser irradiation and relate it to the observed alteration in the spectra. Methods. We used olivine (OL) and pyroxene (OPX) pellets irradiated by 5 keV H+ ions and individual femtosecond laser pulses and measured their visible (VIS) and near-infrared (NIR) spectra. We observed the pellets with scanning and transmission electron microscopy. We studied structural, mineralogical, and chemical modifications in the samples. Finally, we connected the material observations to changes in the reflectance spectra. Results. In both minerals, H+ irradiation induces partially amorphous sub-surface layers containing small vesicles. In OL pellets, these vesicles are more tightly packed than in OPX ones. Any related spectral change is mainly in the VIS spectral slope. Changes due to laser irradiation are mostly dependent on the material's melting temperature. Of all the samples, only the laser-irradiated OL contains nanophase Fe particles, which induce detectable spectral slope change throughout the measured spectral range. Our results suggest that spectral changes at VIS-NIR wavelengths are mainly dependent on the thickness of (partially) amorphous sub-surface layers. Furthermore, amorphisation smooths micro-roughness, increasing the contribution of volume scattering and absorption over surface scattering. Conclusions. Soon after exposure to the space environment, the appearance of partially amorphous sub-surface layers results in rapid changes in the VIS spectral slope. In later stages (onset of micrometeoroid bombardment), we expect an emergence of nanoparticles to also mildly affect the NIR spectral slope. An increase in the dimensions of amorphous layers and vesicles in the more space-weathered material will only cause band-depth variation and darkening.Peer reviewe

Spin vector and shape of (6070) Rheinland and their implications

Main belt asteroids (6070) Rheinland and (54827) 2001NQ8 belong to a small population of couples of bodies which reside on very similar heliocentric orbits. Vokrouhlicky & Nesvorny (2008, AJ 136, 280) promoted a term "asteroid pairs", pointing out their common origin within the past tens to hundreds of ky. Previous attempts to reconstruct the initial configuration of Rheinland and 2001NQ8 at the time of their separation have led to the prediction that Rheinland's rotation should be retrograde. Here we report extensive photometric observations of this asteroid and use the lightcurve inversion technique to directly determine its rotation state and shape. We confirm the retrograde sense of rotation of Rheinland, with obliquity value constrained to be >= 140 deg. The ecliptic longitude of the pole position is not well constrained as yet. The asymmetric behavior of Rheinland's lightcurve reflects a sharp, near-planar edge in our convex shape representation of this asteroid. Our calibrated observations in the red filter also allow us to determine $H_R = 13.68\pm 0.05$ and $G = 0.31\pm 0.05$ values of the H-G system. With the characteristic color index $V-R = 0.49\pm 0.05$ for the S-type asteroids, we thus obtain $H = 14.17\pm 0.07$ for the absolute magnitude of (6070) Rheinland. This a significantly larger value than previously obtained from analysis of the astrometric survey observations. We next use the obliquity constraint for Rheinland to eliminate some degree of uncertainty in the past propagation of its orbit. This is because the sign of the past secular change of its semimajor axis due to the Yarkovsky effect is now constrained. Determination of the rotation state of the secondary component, asteroid (54827) 2001NQ8, is the key element in further constraining the age of the pair and its formation process.Comment: Published in AJ, 28 pages, 4 figures, 2 table

Asteroid (16) Psyche’s primordial shape: A possible Jacobi ellipsoid

Context. Asteroid (16) Psyche is the largest M-type asteroid in the main belt and the target of the NASA Psyche mission. It is also the only asteroid of this size (D >  200 km) known to be metal rich. Although various hypotheses have been proposed to explain the rather unique physical properties of this asteroid, a perfect understanding of its formation and bulk composition is still missing. Aims. We aim to refine the shape and bulk density of (16) Psyche and to perform a thorough analysis of its shape to better constrain possible formation scenarios and the structure of its interior. Methods. We obtained disk-resolved VLT/SPHERE/ZIMPOL images acquired within our ESO large program (ID 199.C-0074), which complement similar data obtained in 2018. Both data sets offer a complete coverage of Psyche’s surface. These images were used to reconstruct the three-dimensional (3D) shape of Psyche with two independent shape modeling algorithms (MPCD and ADAM). A shape analysis was subsequently performed, including a comparison with equilibrium figures and the identification of mass deficit regions. Results. Our 3D shape along with existing mass estimates imply a density of 4.20  ±  0.60 g cm−3, which is so far the highest for a solar system object following the four telluric planets. Furthermore, the shape of Psyche presents small deviations from an ellipsoid, that is, prominently three large depressions along its equator. The flatness and density of Psyche are compatible with a formation at hydrostatic equilibrium as a Jacobi ellipsoid with a shorter rotation period of ∼3h. Later impacts may have slowed down Psyche’s rotation, which is currently ∼4.2 h, while also creating the imaged depressions. Conclusions. Our results open the possibility that Psyche acquired its primordial shape either after a giant impact while its interior was already frozen or while its interior was still molten owing to the decay of the short-lived radionuclide 26Al.Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 199.C-0074 (principal investigator: P. Vernazza). P. Vernazza, A. Drouard, M. Ferrais and B. Carry were supported by CNRS/INSU/PNP. J.H. and J.D. were supported by grant 18-09470S of the Czech Science Foundation and by the Charles University Research Programme no. UNCE/SCI/023. E.J. is F.R.S.-FNRS Senior Research Associate. The work of TSR was carried out through grant APOSTD/2019/046 by Generalitat Valenciana (Spain). This work was supported by the MINECO (Spanish Ministry of Economy) through grant RTI2018-095076-B-C21 (MINECO/FEDER, UE)

Binary asteroid (31) Euphrosyne: ice-rich and nearly spherical

Aims. Asteroid (31) Euphrosyne is one of the biggest objects in the asteroid main belt and it is also the largest member of its namesake family. The Euphrosyne family occupies a highly inclined region in the outer main belt and contains a remarkably large number of members, which is interpreted as an outcome of a disruptive cratering event. Methods. The goals of this adaptive-optics imaging study are threefold: to characterize the shape of Euphrosyne, to constrain its density, and to search for the large craters that may be associated with the family formation event. Results. We obtained disk-resolved images of Euphrosyne using SPHERE/ZIMPOL at the ESO 8.2 m VLT as part of our large program (ID: 199.C-0074, PI: Vernazza). We reconstructed its 3D shape via the ADAM shape modeling algorithm based on the SPHERE images and the available light curves of this asteroid. We analyzed the dynamics of the satellite with the Genoid meta-heuristic algorithm. Finally, we studied the shape of Euphrosyne using hydrostatic equilibrium models. Conclusions. Our SPHERE observations show that Euphrosyne has a nearly spherical shape with the sphericity index of 0.9888 and its surface lacks large impact craters. Euphrosyne’s diameter is 268 ± 6 km, making it one of the top ten largest main belt asteroids. We detected a satellite of Euphrosyne – S/2019 (31) 1 – that is about 4 km across, on a circular orbit. The mass determined from the orbit of the satellite together with the volume computed from the shape model imply a density of 1665 ± 242 kg m−3, suggesting that Euphrosyne probably contains a large fraction of water ice in its interior. We find that the spherical shape of Euphrosyne is a result of the reaccumulation process following the impact, as in the case of (10) Hygiea. However, our shape analysis reveals that, contrary to Hygiea, the axis ratios of Euphrosyne significantly differ from those suggested by fluid hydrostatic equilibrium following reaccumulation.This work has been supported by the Czech Science Foundation through grant 18-09470S (J. Hanuš, O. Chrenko, P. Ševeček) and by the Charles University Research program No. UNCE/SCI/023. M.B. was supported by the Czech Science Foundation grant 18-04514J. Computational resources were supplied by the Ministry of Education, Youth and Sports of the Czech Republic under the projects CESNET (LM2015042) and IT4Innovations National Supercomputing Centre (LM2015070). P. Vernazza, A. Drouard, M. Ferrais and B. Carry were supported by CNRS/INSU/PNP. M.M. was supported by the National Aeronautics and Space Administration under grant No. 80NSSC18K0849 issued through the Planetary Astronomy Program. The work of TSR was carried out through grant APOSTD/2019/046 by Generalitat Valenciana (Spain). This work was supported by the MINECO (Spanish Ministry of Economy) through grant RTI2018-095076-B-C21 (MINECO/FEDER, UE). The research leading to these results has received funding from the ARC grant for Concerted Research Actions, financed by the Wallonia-Brussels Federation. TRAPPIST is a project funded by the Belgian Fonds (National) de la Recherche Scientifique (F.R.S.-FNRS) under grant FRFC 2.5.594.09.F. TRAPPIST-North is a project funded by the Université de Liège, and performed in collaboration with Cadi Ayyad University of Marrakesh. E. Jehin is a FNRS Senior Research Associate